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Abstract. The aim of the Yakutsk array enhancement project 
is to create an instrument to study the highest-energy galac- 
tic cosmic rays (CRs) - their sources, energy spectrum, and 
mass composition. Additionally, there will be unique ca- 
pabilities for investigations in the transition region between 
galactic and extragalactic components of CRs. 

Using the well-developed imaging atmospheric 
Cherenkov telescope technique adapted to the energy 
region E > 10 15 eV, we plan to measure the longitudinal 
structure parameters of the shower, e.g., angular and tem- 
poral distributions of the Cherenkov signal related to X max 
and the mass composition of CRs. 

The main advantages of the Yakutsk array, such as its 
multi-component measurements of extensive air showers, 
and model-independent CR energy estimation based on 
Cherenkov light measurements, will be inherited by the in- 
strument to be created. 



1 Introduction 

The cosmic ray (CR) astrophysics behind extensive air 
shower (EAS) measurements in the energy range E > 10 14 
eV pose several problems of great interest that stimulate 
further in-depth research. In addition to general objectives 
such as CR energy spectrum, mass composition, and sources, 
some specific questions should be addressed in the immedi- 
ate future. For example: What is the maximum energy of 
supernova remnant (SNR) accelerators? Is the 'knee' of CR 
spectrum due to the diffusion of particles in magnetic fields 
or to the upper limit of galactic sources? Where is the tran- 
sition region between galactic and extragalactic components 
of CRs? 

Some theoretical predictions may soon be verified exper- 
imentally. In particular, it is shown in nonlinear kinetic the- 




Fig. 1. Plan of the Yakutsk array. Captions: open circles = soft 
component detectors; filled circles = Cherenkov light detectors; 
squares = muon detectors. Three stations at the bottom (indicated 
by crosses) were destroyed by a flood in 2010. 



ory of CR acceleration in SNRs that t he maximum energy 
of nu clei is proportional to their charge ( iBerezhko and VoUa 



2007). Thus, the knee in the energy spectrum is the result of 
the contribution of progressively heavier species, and galac- 
tic CRs are dominated by iron-grouj ) nuclei at E ~ 10 17 eV. 



An alt ernative model is proposed bv lErlvkin and W olfendale 
(2001) in which a recent nearby SNR accelerates CRs re- 
sulting in rather different structure in the energy spectrum, 
detectable as well. 

As an initial step in enhancement of CR detecting possi- 
bilities, the Yakutsk array was recently equipped with a pin- 
hole camera with a row of photomultiplier tubes (PMTs) in 
the rear as a prototype detector of the angular distribution 
of Cherenkov light i nduced by showers in the atmosphere 
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dGaripov et al.L 1200 II) . A conclusion from the measurements 
was that angular and temporal characteristics of the signal 
can be used to investigate longitudinal shower profile in the 
atmosphere. 

Developing these efforts further, we plan to use new device 
consisting of multi-anode PMT at the focus of a spherical 
mirror as a Cherenkov light detector embedded in the array. 
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2 Present status of the Yakutsk array, and the next-step 
astrophysical goals 

At present, the Yakutsk array is measuring the soft compo- 
nent of EAS with 58 stations, muons with 6 underground 
detectors, and Cherenkov light with 48 PMTs. All detectors 
are irregularly distributed within a 10 km 2 array area (Fig.[TJ; 
the target energy range of investigations is 10 15 eV to 10 20 
eV. More technical det a ils and physica l resul ts are given in 



Dvakonov et al .l (119911) ; Egorova et al. (2004J); llvanov et al 
d2009l) : llvanov|j2010l) and lthe website! (120061) . 



Our next task is to modernize the array in order to ob- 
tain a precise instrument capable of measuring the highest- 
energy galactic CRs - their sources, energy spectrum, and 
mass composition. Related to this, we also aim to study a 
transition region between galactic and extragalactic compo- 
nents of CRs where some irregularities in spectrum and com- 
position may be revealed. 

A crucial clue in this approach is the accurate determi- 
nation of the mass composition of CRs, which is a weak 
point of existing giant EAS arrays. In this context we plan 
to adapt the we ll-known atmo spheric Cherenkov telescope 
technique (e.g., iHintonl |2009) to measure the angular and 
temporal structure of the signal connected to EAS longitudi- 
nal profiles above E — 10 15 eV. 

The idea is not to concentrate on discriminating the 
gamma-ray initiated showers from the hadronic showers; in- 
stead, all showers from different primary particles will be 
detected by the wide field of view (FoV) telescopes in co- 
incidence with the array detectors. The angular and tempo- 
ral parameters of showers measured will then be analyzed to 
identify EAS primary particles. Longitudinal shower profile 
measurements can also be used to estimate some hadronic 
interaction parameters (IKnurenko et al.l 120081) . 

Experimental arguments in elucidating the origin of the 
knee and ankle in CR spectrum will significantly strengthen 
due to the estimation of mass composition in the energy 
range relevant. Existing models (e.g., listed in Introduction) 
are very different in composition expected around the knee 
and ankle, so the estimation of the average mass of the pri- 
maries in addition to the improved measurement of the sharp- 
ness of the knee/ankle should allow us to discriminate be- 
tween models. 

Another important step in the adaptation is the essen- 
tial reduction in the size and cost of the telescope. This 
is possible through increased threshold energy because the 
total number of Cherenkov photons emitted by EAS rela- 
tivistic electrons is proportional to the shower energy. If 
we compare one of the High-Energy Stereoscopic System 
(HESS) telescopes (diameter D = 13 m, E thr = 10 11 eV 
/www.mpi-hd.mpg.de/hfm/HESS/HESS.shtml/) with that re- 
duced to D = 13 cm, then the number of Cherenkov photons 
detected from EAS is comparable at E > 10 15 eV. Hence, 
in our energy range, the set of Cherenkov telescopes with 
D = 13 cm is approximately equivalent to HESS functioning 



in the energy range E > 10 11 eV; this is, of course, except- 
ing the event rate: CR intensity ratio in two energy intervals 
is J(E > 10 15 ) /J(E > 10 11 ) = 10~ 4k = 1.6 x 10" 7 , where 
K — 1.7 is the integral spectrum index. 



3 Modeling the measurement of EAS longitudinal 
structure parameters with Cherenkov telescopes 

It is well known that the angular and temporal struc- 
ture of Cherenkov light emitted by EAS can be used 
to infer the longitudinal development parameters of a 



shower dV . Zatsepin, 1964; Fomin and Khristiansenl 11971 
iGrigoriev et. all 1 1979b . The angular distri bution of a 
Chere nkov signal from EAS was calculated by lv. Zatseplnl 
(1964), where it was assumed that the angular distribution 
of emitt ed photons is determined by th at of electrons in the 
shower. iFomin and Khristiansen ( 197 lb then proposed using 
the pulse shape of the Cherenkov signal, namely the pulse 
width, to indicate the shower maximum position X max in 
the atmosphere. 

Experimental measurements of the Cherenkov signal pulse 
shape were performed in Yakutsk in 1973-1979. The 
pulse width was measured in the core distance interval 
(100,800) m for t he sho wer size 1.9 x 10 7 < N < 1.3 x 10 8 
( Dvakonov et all 1991). The results were used not only 
to estimate X max , but attempts were also made to eval- 
uate the EAS cascade curve at energie s arou nd 10 17 eV 
dGrigorjev et. all 1 19791: IKnurenko et~all l200ll) . Unfortu- 
nately, the angular and temporal resolutions of optical in- 
struments and data-acquisition systems at that time were not 
sufficient to obtain the accuracy required to distinguish EAS 
primaries and/or nuclear interaction models. 

To verify the ability of the Cherenkov telescope to discrim- 
inate angular and temporal profiles of signals from showers 
initiated by primary nuclei and photons at energies above 
10 15 eV, we mode l ed the process. The same assumptions as 



m lDvakonov et all d 199 lb are used to describe the Cherenkov 
light emission in atmosphere by the vertical EAS: The energy 
spectrum of electrons is taken at s = 1 (shower maximum); 
the spatial distribution of electrons is neglected; angular dis- 
tributions of electrons and Cherenkov photons are Gaussian; 
and the time of flight to detector for photons is given by 
t — — \Jh 2 + R 2 — h/c, where n is the refraction coefficient 
of air, h is height, and R is axis distance. 

In this crude approximation the number of photons is 
given by the integral of the EAS cascade curve multi- 
plied by the function giving the contribution of an elec- 
tron to the Chere nkov light flux at a giv en R and depth 
in the atmosphere (Dvakono v et al. , 199 lb . As the cascade 
curve, we used t he gamma-distribu t ion ap proximation of 
the HiRes results dAbu-Zavvad et all 1200 lb . In the energy 
range (10 15 ,10 17 ) eV beyond the reach of HiRes, we used 
neXus 2 model simulation results for showers initiated by 
P,Fe dKnapp et all l200lb . A difference in X max (E) of 
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Fig. 2. Angular spread of EAS Cherenkov light in the detector 
at R — 800 m from the axis. Vertical showers (E = 10 ° eV) are 
shown. The areas within dotted curves represent the intrinsic fluc- 
tuations of the showers. 



Fig. 3. Temporal structure of the Cherenkov signal from the vertical 
EAS, E = 10 15 eV, R = 800 m, induced by the primary particles 
indicated. Symbols are the same as in the previous Figure. 



showers from proton and iron nuclei in this model was also 
employed as the input parameter for the cascade curve from 
HiRes. For the 7-initiated showers, the Greisen formula was 
used. 

In Fig. 12 the angular distribution of Cherenkov light in 
the detector at R — 800 m is shown. The EAS primaries are 
P,Fe, and 7 with E = 10 15 eV. It appears from this that the 
angular resolution ~ 1° of the Cherenkov telescope is suffi- 
cient to distinguish the primaries (E > 10 15 eV), at least the 
7-quanta from the nuclei, based on the angular distribution 
of the signal alone. 

Another parameter that has been suggested by the pio- 
neers is the temporal characteristic of the Cherenkov sig- 
nal from EAS. In Fig. [3] the results of our simulations are 
shown. In this case the approximation used is too rough, so 
that the half width of the signal is overestimated by ~ 40% 
in comparison with the actual value measured in Yakutsk 
dDvakonov et all 1 1991 ) and with more precise simulations 
(e.g.. iKalmykov et all 1975 ). Nevertheless, assuming that 
the overestimation is independent of the primary particle type 
(diminishing the widths of the signals to 0.7 in all three 
cases), we can conclude that the time resolution St < 10 ns is 
sufficient to distinguish the primaries at E > 10 15 eV. 

Shower fluctuations characterized by oxmax — 67 '/22/47 
g/cm 2 ,forP/Fe/7 primaries respectively, are shown by the 
dotted curves in Figs [2] [3] Conclusions about the possibility 
to distinguish primaries do not weaken because of fluctua- 
tions less than the differences caused by the primaries. More 
reliable results will be achieved by combining both the angu- 
lar and the temporal measurements in one experiment. 



4 Preliminary Cherenkov telescope design, and en- 
hancement of the Yakutsk array 

Recent developments in the PMT industry have made it pos- 
sible to dramatically reduce the size of Cherenkov telescopes. 



For example, the imaging camera for the HESS telescope 
consists of 960 PMTs each with a photocathode area of ~ 16 
cm 2 (pixel size). Proportionally reducing the area of mir- 
ror and camera 10 4 times, we should have a pixel size of 
~ 0.4 x 0.4 mm 2 for the same angular resolution. 

A newly developed semiconductor light-sensor, the so- 
called Geiger-mode avalanche photodiode (G-APD, or 'sil- 
icon PMT'), typically has a sensitive area of ~ 1 x 1 mm 2 
covered by 100 to 10 4 cells. This is larger than we need, but 
on the other hand, the mirror diameter can be greater than 13 
cm. Application of this type of light- sensor is currently bei ng 
widely discussed - see for example iTeshima et all ( 12007b - 
and there are examples of cam eras designed w i th app ropriate 
data-acquisition systems (e.g., Anterhub et al. . 20091) . 

Another type of light-sensor applicable to the Cherenkov 
mini-telescope is position-sensitive PMT. In this case the flat 
multi-anode is segmented innxn pixels (n = 2, ..,16 for 
Hamamatsu products; pixel area > 2 x 2 mm 2 ). 

We selected Hamamatsu R2486 PMT with a 16 x 16 
crossed wire anode for the prototype detector. The PMT is 
used as an imaging camera placed near the focus of a spheri- 
cal mirror. Because the radius of the R2486 photocathode is 
3 cm, the mirror diameter was chosen to be D = 26 cm, and 
the radius of curvature R = 22.5 cm. The mirror parameters 
were optimized to provide as wide a FoV as possible where 
the image distortions are less than the given pixel size. 

Angular resolution of the prototype telescope is deter- 
mined by the pixel size d = 3.8 mm and focal distance 6a — 
2d/i?~2°. The width of FoV was estimated as (-14°, 14°) 
in regard to the optical axis. 

The result of modeling in the plane intersecting the mirror 
optical axis is shown in Fig. [4] The light intensity emitted by 
the distant point sources with slant angles —14° < a < 14° 
is shown as a function of the axis distance on the camera 
surface. The root-mean-square radius of the light spot av- 
eraged over camera surface is 0.64 of the pixel size. Thus, 
the position of intensity maximum is an explicit function of 
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Fig. 4. Scanning angular distribution of light from a distant point 
source at the slant angle a with the position-sensitive PMT in the 
focus of the spherical mirror. 



incident angle. This demonstrates the ability of our proto- 
type Cherenkov telescope to convert the angular distribution 
of light sources to the spatial distribution of signals on the 
imaging camera surface with adequate accuracy. 

The effective acceptance area of our prototype telescope 
is ~ 200 cm 2 taking into account the device geometry, shad- 
owing of the mirror, etc. This is comparable to the area of 
FEU-49 PMT (176 cm 2 ) in use in Cherenkov detectors of 
the Yakutsk array. Thus, the real EAS event rate previously 
detected with the Cherenkov sub-array can be used as an esti- 
mator of the expected number of showers. In the arrangement 
of Cherenkov light detectors shown in Fig.Q] the number of 
showers detected per hour at£'>1.2xl0 15 eV is 62, and 
4 eve nts were observed at E > 6 x 10 16 eV dlvanov et all 
20031) . 

If we place the set of Cherenkov telescopes in a congru- 
ent hierarchical triangle structure as in Fig. Q] (filled circles), 
then the expected number of EASs detected per year (as- 
suming 500 hours/year is suitable for Cherenkov light mea- 
surements) is ~ 31000 at E > 1.2 X 10 15 eV, and - 2000 at 
£>6x 10 16 eV. 



5 Conclusions 

We have demonstrated the possibility of investigating EAS 
longitudinal development in the atmosphere through measur- 
ing the angular and temporal distributions of a Cherenkov 
signal with a telescope embedded in a surface EAS array en- 
vironment. The enhancement of the Yakutsk array by a set of 
Cherenkov telescopes adapted to the energy range E > 10 15 
eV is aimed at significantly improving the array's ability to 
study the mass composition of CRs. The main advantage 
of this combination consists in the integration of the multi- 
component measurement of EAS and model-independent CR 



energy estimation by the surface array with excellent angular 
and temporal resolution of the Cherenkov telescopes. 
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